The effect of magnetic impurities in a two-band superconductor: 
A point-contact study of Mn-substituted MgB 2 single crystals 
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We present the first results of directional point-contact measurements in Mgi_ x Mn x B2 single 
crystals, with x up to 0.015 and bulk T c down to 13.3 K. The order parameters A CT and A,r were 
obtained by fitting the conductance curves with the two-band Blonder- Tinkham-Klapwijk (BTK) 
model. Both A T and A CT decrease with the critical temperature of the junctions T^, but remain 
clearly distinct up to the highest Mn content. Once analyzed within the Eliashberg theory, the results 
indicate that spin-flip scattering is dominant in the a band, as also confirmed by first-principle band 
structure calculations. 



The two-band character of superconductivity in MgB2 
Q, has been almost completely understood and ex- 
plained by now, but the effects of disorder and chemical 
doping arc still in need of some experimental clarification. 
Several substitutions have been tried, either in the Mg or 
in the B sites, but only few of them have been successful 
and none has been able to enhance the critical temper- 
ature of the compound. However, their study has proven 
useful to clarify the role of the different scattering chan- 
nels and to try to control them selectively. Substitutions 
with magnetic impurities (Mn, Fe) represent a class of its 
own because of the spin-flip pair-breaking scattering that 
is expected to dramatically suppress superconductivity - 
even though the way it does it in a two-band supercon- 
ductor has never been studied experimentally Up to 
now, only two successful Mn substitution in the Mg site 
have been reported, one in polycrystalline samples [4| and 
one in single crystals 0. Here, we present the first re- 
sults of point-contact spectroscopy (PCS) measurements 
in Mn-substituted MgB2 single crystals grown at ETHZ, 
which allowed us to study the effects of magnetic im- 
purities on the order parameters (OPs) of a two-band 
superconductor. The amplitudes of the OPs, A a and 
A,,-, were determined as a function of the critical tem- 
perature. The resulting trend can be explained within 
the Eliashberg theory as being mainly due to a doping- 
induced increase in the pair-breaking scattering within 
the a bands, with minor contributions from the 7r-7r or 
the (j-7r channels. This result, apparently in contrast 
with the Mn position in the lattice, is however confirmed 
by first-principle calculations of the local effect of a Mn 
impurity on the bandstructure of MgB2. 

The high-quality Mgi_ x Mn x B2 single crystals used for 
our measurements were grown by using the same high- 
pressure, cubic-anvil technique set up for pure MgB2, 
and by replacing part of the Mg precursor with metallic 
Mn |5j. The Mn content x of each crystal was measured 
by EDX through a careful evaluation of the Mn/Mg ra- 
tio; the crystals are single-phase and homogeneous within 



±Sx (5x=0.0010). The bulk critical temperature T c was 
determined by DC magnetization measurements. As re- 
ported elsewhere the dependence of the lattice con- 
stants on the Mn content indicates that Mn replaces 
Mg in the lattice. The magnetic-field dependence of 
the Curie part C* of the magnetic moment M clearly 
indicates that Mn ions are divalent (i.e. Mn 2+ ) and 
in the low-spin state (S — 1/2), as also confirmed by 
our first-principle band-structure calculations. The crys- 
tals we used were carefully selected among those with 
the sharpest superconducting transitions and best struc- 
tural properties, so that secondary phases or other im- 
purities can be excluded. They had different Mn con- 
tents x between 0.0037 and 0.0150, corresponding to bulk 
critical temperatures T c between 33.9 and 13.3 K, with 
AT c (10-90%) increasing with doping from 0.65 K to 5.4 
K. The point contacts were made by putting a small drop 
of silverpaint (0 ~ 50^m) on the flat side surface of the 
crystal |6y so as to inject the current mainly parallel to 
the ab planes, which in pure MgB2 is the best configura- 
tion for the observation of both the a- and 7r-band gaps 
0, El- In most cases, we studied the temperature de- 
pendence and the magnetic-field dependence of the con- 
ductance curves (dl/dV vs. V), so as to determine the 
critical temperature of the junction (i.e., the "Andreev 
critical temperature", T<r), and to understand whether 
one or two OPs were present 0, @- The conductance 
curves were then normalized to the normal state and fit- 
ted with the two-band BTK model || , as described else- 
where 0. The fitting function contains the OPs A CTi7r , 
the coefficients Z a ^ (related to the potential barrier and 
to the the Fermi velocity mismatch at the interface) and 
the lifetime broadening coefficients T „^ as adjustable pa- 
rameters, plus the weight w n of the 7r-band contribution 
to the conductance. For reliable estimates of the OP 
amplitudes, whenever possible, we selected contacts with 
rather high resistance (i?N ^ 30 O) , and with no dips 
in the conductance curves, so as to fulfill the requirements 
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FIG. 1: Normalized conductance curves (symbols) of different 
afr-plane junctions in Mgi_ x Mn x B2 single crystals at T—4.2 
K. The curves are labeled with the corresponding T C A values. 
The Mn content and the normal-state junction resistance are 
shown as well. Lines: best-fit curves given by the two-band 
BTK model. The 7r-band weight is also indicated. In 
the absence of specific theoretical predictions, we took 0.66 < 
W-,, < 0.8 as for afc-plane contacts in pure MgB2 QUI. 

for ballistic conduction |To| . 

Fig. ^ reports some normalized conductance curves 
(symbols) in crystals with different Mn content. From 
now on, we will label the curves with the corresponding 
value of T A instead of the Mn content or the bulk T c , 
since PCS is a local, surface-sensitive probe. T C A can be 
smaller than T c if the bulk transition is broad (as for the 
most-doped crystals) or if proximity effect occurs at the 
S/N interface (when £ « a, being a the true contact size). 

Lines in Fig. ^ represent the two-band BTK fit of 
experimental curves. In conventional superconductors, 
spin- flip scattering makes the superconducting gap be- 
come ill-defined and, for strong scattering, bands of 
states within the original energy gap are formed |l2l Il3| . 
Strong-coupling calculations for MgB 2 01 nave shown 
that magnetic impurities can give rise to sub-gap states 
in both the partial a- and rr-band DOS, that are not 
taken into account by the BTK model. However, the 
sub-gap features in the DOS are by far smaller than the 
peaks connected to the OPs A T and A CT ^ij. At any 
finite temperature they are further smeared out, so that 
they would be very difficult to observe experimentally. 
One can reasonably expect that the same happens in 
the Andreev-reflection curves. As a matter of fact, all 
the experimental conductance curves are rather broad- 
ened and have smaller amplitude (see Fig^l than in pure 
MgB 2 :6|, so that even if sub-gap structures exist, they 
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FIG. 2: Experimental conductance curves at T=4.2 K of two 
contacts having T C A =3I K (a) and T C A =I8 K (b), in the pres- 
ence of different magnetic fields B || c (circles) and the rele- 
vant two-band BTK fit (lines). Vertical dotted lines indicate 
the position of the conductance peaks. The values of the fit- 
ting parameters in zero field are: A^-=2.5 meV, Ao-=5.2 meV; 
T^=2.3 meV, r CT =6.0 meV; ^=0.55, Z CT =0.82, to^=0.7 for 
panel (a); A^=1.56 meV, A CT =2.8 meV; r^=2.0I meV, 
F CT =3.8 meV; ^=0.43, Z CT =0.50, 1^=0.7 for panel (b). The 
two-band BTK fit works well up to the critical field B C 2 (— 
3.5 T in (a) and ~ 2.5 T in (b)). 

are practically unobservable. Hence, the BTK model can 
be used as a reasonable, first-order approximation to a 
more specific model for Andreev reflection in the pres- 
ence of magnetic impurities on the superconducting side 
of the contact, which is lacking at the present moment. 

Even if the fit of the conductance curves shown in Fig.Q] 
indicates the existence of two OPs in the Mgi_ x Mn x B2 
system, it is clear that some curves - especially at the 
highest doping levels - show little or no structures associ- 
ated to the larger OP. To check in a more convincing way 
if two-band superconductivity persists up to the highest 
Mn content, we studied the magnetic-field dependence of 
the dl/dT^ curves in the whole range of T C A . 

Fig. Ufa) shows the magnetic-field dependence of the 
conductance curve of a contact with T C A =31 K (circles) 
with the relevant two-band BTK fit (lines). The zero- 
field curve shows both the peaks corresponding to A,,- 
and the smooth shoulders related to A CT . On apply- 
ing the magnetic field, the small-gap features reduce in 
amplitude, progressively unveiling the underlying large- 
gap features. An outward shift of the peaks is observed 
at some B = B*, when the tr-band structures become 
dominant and determine the shape of the conductance 
curve. Unlike in pure MgB2, here B* is intense enough 
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FIG. 3: Symbols: zero-field, low-temperature (T=4.2 K) con- 
ductance curve of a contact with T<* = 9.6 K made on the 
most-doped crystals [x = 0.015). Solid line: best fit given 
by the two-band BTK model. Dashed line: best fit given by 
the standard (single-band) BTK model. Arrows indicate the 
smooth structures related to the two order parameters, which 
are magnified in the inset for clarity. 

(if compared to the critical field, strongly suppressed by 
Mn impurities) to partly depress A CT , while the 7r-band 
partial conductance seems not to vanish at B* . This 
makes it impossible to separate the partial a- and 7r- 
band conductances as we did in pure MgB 2 0, , but 
the shift of the conductance peaks is anyway a clear proof 
of the existence of two OPs. The same happens in the 
magnetic-field dependence of contacts with much lower 
critical temperature. For example, Fig. Gib) reports the 
field dependence of the normalized conductance curve of 
a contact with T C A =18 K. Here the zero-field conductance 
does not show clear structures related to A CT but an out- 
ward shift of the conductance maxima is observed any- 
way at B = B* ~ 1 T. For higher doping levels, i.e. for 
T C A < 17 K, the low-temperature, zero-field conductance 
curves can be very well fitted with the two-band BTK 
model (see the bottom curve in FigCJ, but the critical 
field is so small and the conductance curves are so broad- 
ened that their magnetic-field dependence is not conclu- 
sive, at least with our experimental resolution. Moreover, 
some conductance curves in this region allow a standard 
(i.e. single-band) BTK fit as well, which gives an OP am- 
plitude A that is the "average" of the values given by the 
two-band fit. However, the single-band fit often fails in 
reproducing both the position of the peaks and the width 
of the Andreev-reflection structures in the conductance, 
as shown in Fig.[3]for a contact with the lowest T A . This 
suggests that two OPs are likely to be present also in the 
most doped samples. 

The complete dependence of the OPs on T C A is reported 
m Fig.H Vertical error bars include both the uncertainty 
on the OP in each fit and the spread of OP values over 
different contacts with the same T C A (within the experi- 
mental uncertainty on T A represented by horizontal error 
bars). For T C A < 1 7 K only the results of the two-band 
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FIG. 4: Symbols: amplitudes of the order parameters as 
a function of the critical temperature of the junctions, T<* . 
Lines: theoretical curves obtained by solving the Eliashberg 
equations in the presence of magnetic scattering in the cr-er 
and tt-tv channels (case (a), solid lines) or in the a-a and <j-tt 
channels (case (b), dashed lines). 

BTK fit are shown, based on the preceding discussion 
and on the regular trend of the OPs for T C A > 18 K. In 
fact, the presence of a single OP in the low-T c A region 
would imply a sudden change in the slope of the curves 
around T C A = 17 K that is not justified by any observed 
discontinuity in the physical properties of the compound 
U . The persistence of two distinct OPs up to the highest 
x value indicates that Mn doping does not significantly 
increase non-spin-flip interband scattering. Actually, this 
is also suggested by the very low Mn concentration (1.5 
% at most) and by the decrease of both the OPs with 
T A (while interband scattering would increase A^[l5|). 
This decrease could arise from changes in the DOSs at 
the Fermi level and in the phonon spectra o?F a ,n{w), 
but these two effects are certainly negligible here, because 
Mn is isovalent with Mg [j| and the Mn content is very 
small. The latter reason also allows assuming the effect 
of ct and 7r intraband non-spin-flip scattering to be much 
smaller than that of pair-breaking scattering in the same 
channels. This leads to conclude that the main possible 
cause of the experimental trend of Fig. 01 is an increase in 
the spin-flip scattering in the ct-ct, tt-tt or ct-7t channels. 

This simple picture can be made quantitative by solv- 
ing the Eliashberg equations (EE) in the presence of ran- 
domly distributed magnetic impurities, treated within 
the Born approximation |l6j . According to the reason- 
ing above, we used the same phonon spectra, DOS val- 
ues and Coulomb pseudopotential as in pure MgB 2 0>@1j 
and neglected all the non-spin-flip scattering rates. We 
thus took as the only adjustable parameters the spin- 
flip scattering rates within the bands (T aa and F^) and 
between bands (r CTlr ). We immediately found that Y aa 
is necessary to fit the T A and OP values and that it 
must be greater than both Y^n and Y a7 , - otherwise A w 
decreases too fast on decreasing T C A . As for r CT7r , this 
agrees with the predictions of Ref. For simplicity, 
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we analyzed separately the two cases: (a) T acr > T^„, 
Tcttt^O and (b) T aa > r CT7r , T^^O. First of all, we fixed 
r 7r7r =r (T(T =r CT7r =0 at T C A =39.4 K (i.e. x=0). Second, wc 
found the values of the parameters that give the exper- 
imental values of A a and A w at T A =18.0 K (cases (a) 
and (b)), by solving the imaginary- axis EE and analyti- 
cally continuing the solution to the real axis. Once deter- 
mined the values of and T a7r in these two points, we 
searched for the simplest r ww (T c A ) and r CTT (T c A ) curves 
connecting them and allowing the fit of the experimental 
values of T A and of the OPs in the whole doping range, 
with no restrictions on T aa . We found out that these 
curves are a parabola (for r wff ) and a straight line (for 
T a7T ). In both cases, T aa follows an almost parabolic 
trend as a function of T A (see FigEJ. The resulting 
theoretical curves that best fit the OPs are reported in 
Fig^as solid lines (case (a)) and dashed lines (case (b)). 
The agreement between experimental data and theoreti- 
cal calculations (especially in case (a)) is striking in the 
whole range of T C A . That our simple - but reasonable 
- model works so well suggests that, up to x = 0.015, 
Mgi_ x Mn x B2 can be treated as a perturbation of the 
unsubstituted compound, with only the addition of a 
intraband magnetic scattering and minor contributions 
from either the n intraband or the interband scattering. 

The intense Mn pair-breaking scattering in the a-a 
channel has been recently predicted as being due to the 
hybridization of the a bands of MgB 2 with the d or- 
bitals of Mn 01 ■ To demonstrate that scattering in this 
channel is much greater than in the n-n and a-u ones, 
we performed preliminary calculations of the electronic 
structure near a Mn substitutional impurity in a 2 x 2 x 4 
MgB 2 superlattice. The results show that around the 
Fermi level (Ep) there is a spin-down d z i band, respon- 
sible for the Mn magnetic moment, very sensitive to the 
details of structural parameters [17| . The a and 7r bands 
behave quite differently near the Mn impurity: while the 
7r electron spectral density is depleted around Ep (due to 
the TT-d interaction), the a bands have a large amplitude. 
Furthermore, the a bands show a sizeable (> 15 meV) 
exchange splitting near Ep, larger than the supercon- 
ducting gap - which is consistent with the experimen- 
tally observed complete suppression of superconductivity 
at about 2% of Mn [5j. On the other hand, the impurity- 
induced mixing of a and 7r states around Mn (providing a 
qualitative indication of interband scattering), is present 
but not very important. These results can explain the 
quick drop of T c , with the persistence of the two distinct 
gaps, and a larger scattering within the a band [lij . 

In conclusion, we have presented the results of the first 
experimental study of the effects of magnetic impurities 
on the order parameters A a and A n of a two-band su- 
perconductor. We have shown that, in Mg!_ x Mn x B 2 , 
A CT and A T decrease regularly with the critical tempera- 
ture but remain clearly distinct down to the lowest T A . 
Within the Eliashberg theory, this is due to an increase in 
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FIG. 5: Spin-flip scattering parameters used to fit the data 
of Fig^l as a function of T C A . Solid symbols: case (a), r CT<T 
(squares) and (circles). Open symbols: case (b), F aa 
(squares) and r CT7r (triangles). Lines are guides to the eye. 

spin-flip scattering in the a bands on increasing the Mn 
content, with possible minor contributions from the <j-tt 
or the 7T-7T channels. This somehow unexpected conclu- 
sion is also supported by first-principle calculations of the 
bandstructure of MgB 2 in the vicinity of a Mn impurity. 
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2004022024 and the INTAS Project N. 01-0617. V.A.S. 
acknowledges support by RFBR (Proj. N. 06-02-16490). 
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